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Previous work on neutrino emission from proto-neutron stars which employed full solutions of the 
Boltzmann equation showed that the average energies of emitted electron neutrinos and antineutri- 
nos are closer to one another than predicted by older, more approximate work. This in turn implied 
that the neutrino driven wind is proton rich during its entire life, precluding r-process nucleosyn- 
thesis and the synthesis of Sr, Y, and Zr. This work relied on charged current neutrino interaction 
rates that are appropriate for a free nucleon gas. Here, it is shown in detail that the inclusion of the 
nucleon potential energies and collisional broadening of the response significantly alters this conclu- 
sion. Iso-vector interactions, which give rise to the nuclear symmetry energy, produce a difference 
between the neutron and proton single-particle energies AU — U„ — U p and alter the kinematics 
of the charged current reactions. In neutron-rich matter, and for a given neutrino/antineutrino 
energy, the rate for v e + n — > e~ + p is enhanced while u e + p — > + e + is suppressed because 
the Q value for these reactions is altered by ±A{7, respectively. In the neutrino decoupling region, 
collisional broadening acts to enhance both v e and i> e cross-sections and RPA corrections decrease 
the v e cross-section and increase the v e cross-section, but mean field shifts have a larger effect. 
Therefore, electron neutrinos decouple at lower temperature than when the nucleons are assumed 
to be free and have lower average energies. The change is large enough to allow for a reasonable 
period of time when the neutrino driven wind is predicted to be neutron rich. It is also shown that 
the electron fraction in the wind is influenced by the nuclear symmetry energy. 

PACS numbers: 26.50.+X, 26.60.-c, 21.65.Mn, 95.85.Ry 



I. INTRODUCTION 

The neutrino opacity of dense matter encountered in 
core-collapse supernova is of paramount importance to 
the explosion mechanism, potential nucleosynthesis, su- 
pernova neutrino detection and to the evolution of the 
compact remnant left behind. Matter degeneracy, strong 
and electromagnetic correlations, and multi-particle exci- 
tations have all been shown to be important, especially at 
supra-nuclear densities [e.g. 1-8]. Supernova and proto- 
neutron star (PNS) simulations that employ some sub- 
set of these improvements to the free gas neutrino in- 
teraction rates have found that these corrections play 
a role in shaping the t emp oral and spectral aspects of 
neutrino emission 0, 19l4ll|. Much is still uncertain, es- 
pecially because of the approximations one must make 
regarding weak interactions with the dense background 
medium. A specific issue of importance is the difference 
between the average energies of electron neutrinos and 
electron antineutrinos. This difference is largely deter- 
mined by the charged current reactions v e + n — > p + e~ 
and v e + p — > n + e + in neutron-rich matter at densities 
p ~ 10 12 - 10 14 g/cm 3 . 

Recently, one of the authors has shown that an ac- 
curate treatment of mean field effects in simulations of 
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PNS cooling changes the predicted electron fraction in 
the neutrino driven wind (NDW) [l2j relative to simu- 
lations which do not account for mean field potentials 
in nuclear matter [l(], EH, [T3 |. This difference has sig- 
nificant consequences for the nucleosynthesis expected in 
the NDW [e.g. Il5l - [r7| and for neutrino oscillations out- 
side the neutrino sphere [l8|, E3 . In this work, we dis- 
cuss generic aspects of strong interactions that lead to 
a large asymmetry in the charged current reaction rates 
for electron neutrinos and antineutrinos. We also demon- 
strate that this difference manifests itself in potentially 
observable effects on neutrino spectra from supernovae 
and that the difference depends on the assumed density 
dependence of the nuclear symmetry energy. The effect of 
multi-particle excitations and correlations (via the RPA) 
on the charged current response are also explored. 



Neutron-rich matter at densities and temperatures 
relevant to the neutrino sphere of a PNS is charac- 
terized by degenerate relativistic electrons and non- 
relativistic partially degenerate neutrons and protons. 
Beta-equilibrium, with net electron neutrino number 
Yv e — is a reasonably good approximation for the mate- 
rial near the neutrino sphere because, by definition, this 
material can efficiently lose net electron neutrino num- 
ber. At these densities, effects due to strong interactions 
modify the equation of state and the beta-equilibrium 
abundances of neutron and protons. Simple models for 
the nuclear equation of state predict that the nucleon 
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potential energy is 



where z = (qo + (12 — Hi)/T and 



U n / P ~ Vis (n n + n p ) ± V iv (n„ 



n p ) . 



(1) 



where V\ s and Vi v are the effective iso-scalar and iso- 
vector potentials. Empirical properties of nuclear matter 
and neutron-rich matter suggest that V; s x no « —50 MeV 
and Viv x no « 20 MeV. The potential energy associated 
with n — > p conversion in the medium is 



AU = U n - U p 



40 x (n " " p) MeV, (2) 
n 



where no = 0.16 nucleons/fm 3 is the number density at 
saturation. It will be shown that AU changes the kine- 
matics of charged current reactions, so that the Q- value 
for the reaction v e + n — > e~ + p is enhanced by AU 
while that for D e + p — > e + + n is reduced by the same 
amount. The effect is similar to the enhancement due to 
the neutron-proton mass difference, but is larger when 
the number density n > no/20. 

In section UH charged current neutrino opacities in 
an interacting medium are discussed. We consider how 
mean fields affect the response of the medium in detail 
and how this depends on the properties of the nuclear 
equation of state. The effects of nuclear correlations and 
multi-particle hole excitations are also discussed. In sec- 
tion IIII1 the effect of variations of the charged current 
reaction rates on the properties of the emitted neutrinos 
is studied. 



II. THE CHARGED CURRENT RESPONSE 

The differential absorption rate for electron neutrinos 
by the process u e + n — > e~ + p is given by 
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where S T (go, q) and S aT (qo,q) are the response functions 
associated with the Fermi and Gamow- Teller operators, 
r+ and OT-f, respectively. The energy transfer to the 
nuclear medium is qo = E v — E e , and the magnitude of 
the momentum transfer to the medium is q 2 — E 2 + E 2 — 
2E u E e cos 8. In a non-interacting Fermi gas, the response 
functions S T (q Q ,q) = S aT (q ,q) = S F (q ,q) are given by 

M<?o,9) = ^2 / rf3 P2<5(<70 + E 2 - EMI - hi (4) 

where the particle labeled 2 is the incoming nucleon and 
the particle labeled 4 is the outgoing nucleon. When the 
dispersion relation for nucleons is given by E(p) = M + 
p 2 /2M - neglecting the neutron- proton mass difference 
for simplicity - the integrals in Eq. 0] can be performed 
to obtain 
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is the free particle-hole polarization function. fi2 and /14 
are the chemical potentials of the incoming and outgoing 
nucleons, M is the nucleon mass, and 
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e m ; n arises from the kinematic restrictions imposed by 
energy-momentum transfer and the energy conserving 
delta function. Physically, e m i n is the minimum energy 
of the nucleon in the initial state that can accept mo- 
mentum q and energy qo- 



A. Frustrated Kinematics 

The differential cross-section of 24 absorption is the 
product of the nucleon response times the available elec- 
tron phase space 



Pe E e (1 - f e (E e )) « E 2 e exp 



T 
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Due to the high electron degeneracy, the lepton phase 
space increases exponentially with the electron energy. 
To completely overcome electron blocking requires E e — 
E„ e — go ~ Me or go ~ — Me when E Vc <C /i e . How- 
ever, the fermi gas response function in Eq. |4] is peaked 
at qo ~ q 2 /2M ss reflecting the fact that nucleons are 
heavy. At large \qo\ ~ q ~ (x e the response is exponen- 
tially suppressed due to kinematic restrictions imposed 
by Eq. [7] which implies only neutrons with energy 
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can participate in the reaction. For conditions in the PNS 
decoupling region, and in the fermi gas approximation, 
the u e reaction proceeds at q w at the expense of large 
electron blocking. Thus effects that can shift strength 
to more negative qo can increase the electron absorption 
rate exponentially. 

It is well known that the neutron-proton mass differ- 
ence AM = M n — M p increases the Q value for this 
reaction and a more general expression for S(qa,q) de- 
rived in [l[ includes this effect. The effect of AM can 
be understood by noting, that at leading order, it only 
changes the argument of the energy delta-function in Eq. 
[4] and is subsumed by the replacements qo — > {qo + AM) 
and 
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This shift changes the location of the peak of the response 
by moving it to the region where E e is larger and con- 
firming that it increases the Q value and the final state 
electron energy by AM = M n — M p . From Eq. [3] we see 
that the rate for v e absorption is increased by roughly a 
factor (1 + AM/E e ) 2 exp (AM/T). By the same token, 
the Q value for the reaction 9 e + p — > e + + n is reduced 
by AM and this acts to reduce the rate. In this case, the 
detailed balance factor [1 — exp (— z)]^ 1 in the response 
function S(qo, q) is the source of exponential suppression 
- simply indicating a paucity of high energy protons in 
the plasma. For small go <C Pe, the detailed balance 
factor is 

-«exp — - — , (11 

1 - exp (-z) \ T J 

where we have used the fact that p, n — p, p = p, e in beta- 
equilibrium. Since qo —¥ (qo — AM) for the v e process, 
AM will suppress this rate exponentially. This is in line 
with the expectation that AM increases the cross-section 
for v e absorption and decreases it for v e absorption. In 
the following we show that the mean field energy shift, 
driven by the nuclear symmetry energy, has a similar 
but substantially larger effect in neutron-rich matter at 
densities p > 10 12 g/cm 3 . 

B. Mean Field Effects 
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FIG. 1: Top Panel: The electron chemical potential (dashed 
lines) and AU — U n — U p (solid lines) are shown as a function 
of density for the two equation of state models (IUFSU: red 
curves and GM3: black curves) in beta-equilibrium for Y v = 
and T = 8 MeV. The grey band shows an approximate 
range of values for inverse spin relaxation time calculated in 
8] and is discussed in connection with collisional broadening. 
Bottom Panel: The equilibrium electron fraction as a function 
of density for the two equations of state shown in the top 
panel. 



Interactions in the medium alter the single particle en- 
ergies, and nuclear mean field theories predict a nucleon 
dispersion relation of the form 

Ei(k) = Vk 2 + M* 2 + Ui = K(k) + Ui , (12) 

where M* is the nucleon effective mass and Ui is the 
mean field energy shift. For neutron-rich conditions, 
the neutron potential energy is larger due to the iso- 
vector nature of the strong interactions. The difference 
AU = U n — U p is directly related to the nuclear symmetry 
energy, which is the difference between the energy per nu- 
cleon in neutron matter and symmetric nuclear matter. 
Ab-intio methods using Quantum Monte Carlo reported 
in [2(| and [2l| , and chiral effective theory calculations of 
neutron matter by [22| suggest that the symmetry energy 
at sub-nuclear density is larger than predicted by many 
mean field models currently employed in supernova and 
neutron star studies (for a review see 23]). To highlight 
the symmetry energies importance, we choose two mod- 
els for the dense matter equation of state: (i) the GM3 
relativistic mean field theory parameter set without hy- 
perons [24j where the symmetry energy is linear at low 
density; and (ii) the IU-FSU parameter set [25| where the 
symmetry energy is non-linear in the density and large 
at sub-nuclear density. 

The electron chemical potential (dashed lines) and 
neutron-proton potential energy difference (solid lines) 
for these two models are shown as a function of density in 



beta-equilibrium in Figured! Here Y v = for all densities 
and a temperature of 8 MeV is assumed. At sub-nuclear 
densities, the IU-FSU AU is always larger than the GM3 
AU value due to the larger sub-nuclear density symme- 
try energy in the former. The electron chemical potential 
as a function of density, as well as the equilibrium elec- 
tron fraction, is shown in Figure [T] for both models. In 
beta-equilibrium, models with a larger symmetry energy 
predict a larger electron fraction for a given temperature 
and density. Therefore, IU-FSU has a larger equilibrium 
p e than GM3 and the reaction v e + n — > e~ + p will ex- 
perience relatively more final state blocking. However, 
as we show below, the inclusion of AU in the reaction 
kinematics is needed for consistency. 

To elucidate the effects of AU we set M* = M and 
note that this assumption can easily be relaxed [l| and 
it does not change the qualitative discussion below. Be- 
cause in current equation of state models the potential, 
Ui, is independent of the momentum, k, this form of the 
dispersion relation results in a free Fermi gas distribu- 
tion function with single particle energies K(k) for nuclc- 
ons of species i, but with an effective chemical potential 
Pi = pi — Ui. This fact was emphasized in Q, and used 
to show that it was unnecessary to explicitly know the 
values of the nucleon potentials for a given nuclear equa- 
tion of state (which are often not easily available from 
widely used nuclear equations of state in the core-collapse 
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supernova community) when calculating the neutral cur- 
rent response of the nuclear medium. Clearly, if both 
Hi and jii are known, then f/, can be easily obtained. 
This implies that for a given temperature, density and 
electron fraction, the neutral current response function 
is unchanged in the presence of mean field effects, as the 
kinematics of the reaction are unaffected by a constant 
offset in the nucleon single particle energies. In contrast, 
the kinematics of the charged current reaction are af- 
fected by the difference between the neutron and proton 
potentials and the charged current response is altered in 
the presence of mean field effects. 

Inspecting the response function in Eq. [4] and the dis- 
persion relation in Eq. [12] it is easily seen that the mean 
field response is 
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This is obtained from the free gas response by the re- 
placements 



Hi ->■ fa = ^ - Ui 

qa ->■ <Zo = qo + U-2 - U 4 



(16) 



and q —> q. Therefore, we see that the potential difference 
AU = ±(U2 — Ua) affects reaction kinematics and cannot 
be subsumed in the redefinition of the chemical potentials 
(to yield the same individual number densities). 

Because AU > e„ and T for neutrino energies of inter- 
est in the decoupling region, it introduces strong asym- 
metry between the electron neutrino and antincutrino 
charged current interactions because the Q value for the 
reaction v e + n — > e~ + p is increased by AU = U n — U p 
and for v e +P e + + n it is reduced by the same amount. 
Since AU < He, this amount of energy is often not enough 
to put the final state electron above the Fermi surface. 
However, it is enough to put the final state electron in a 
relatively less blocked portion of phase space resulting in 
an exponential enhancement of the cross-section for v e . 
This is shown in Figure [2j where the differential cross- 
section integrated over angle for charged current absorp- 
tion is plotted as a function of the final lepton energy. 
The neutrino energy is set to 12 MeV and the conditions 
of the medium are T — 8 MeV, and ns = 0.02 fm -3 and 
Y e = 0.027. The peak of the differential cross-section is 
shifted by about AU up (down) in e e - (e e +) for elec- 
tron (anti-)neutrino capture. This shift significantly in- 
creases the available phase space for the final state elec- 
tron in v e + n — > e~ + p. The (arbitrarily scaled) phase 
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FIG. 2: Angle integrated differential cross sections for a 
12 MeV neutrino. The solid lines correspond to the reaction 
v e + n — s> e~ +p and the dashed lines correspond to v e +p ^ 
e + + n. The black lines are calculations in which mean field 
effects have been included, while the red lines are calculations 
in which the mean field effects have been ignored. The green 
dotted line corresponds to the available electron phase space, 
arbitrarily scaled. The assumed background conditions are 
T = 8 MeV, and n B = 0.02 fm" 3 . The electron fraction is 
0.027, which corresponds to beta equilibrium for the given 
temperature, density, and the assumed nuclear interactions. 
The nucleon potential difference is U„ — U p = AU = 9 MeV. 
All cross-sections are for the same baryon density and electron 
fraction (i.e. all assume the same p, for the neutrons and 
protons) . 



space factor p e E e (l — f e ) is also plotted and the peak of 
l/Vda/dqo approximately follows this relation. As was 
argued in section lITAl the rate of D e +p — > e+ + n should 
also be approximately proportional to this phase space 
factor and be exponentially suppressed. This is seen in 
the Figure H 

In Figure [31 the inverse mean free path (A -1 = cr/V) 
is shown as a function of neutrino energy for the same 
conditions considered in Figure [2] At low energies the 
electron neutrino mean free path is reduced when mean 
fields are correctly incorporated, but at larger neutrino 
energies the presence of mean fields becomes less impor- 
tant and the mean free paths with and without mean 
fields asymptote to each other. The electron antincutrino 
mean free path is reduced relative to the free gas result 
and the presence of a threshold at the potential difference 
AU is evident in the mean field calculation. The effective 
bremsstrahlung mean free path is also plotted. This is 
calculated assuming the secondary neutrinos are in ther- 
mal equilibrium with the background, which is a good 
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given number density and temperature by inverting the 
free Fermi gas relation, then the response is assumed to 
be the free gas response but with the effective chemi- 
cal potentials in place of the actual chemical potentials. 
This prescription is incorrect because while it accounts 
for the location of the Fermi surface of the nucleons it fails 
to account for the presence of a potential energy differ- 
ence between incoming and outgoing nucleon states. This 
amounts to assuming /j, — »■ jl, so that in Eq. [4] qo — > qa 
and the response becomes the non-interacting response 
for the given density and electron fraction. When the 
potential energies of the incoming and outgoing nucleons 
states are equal, as in symmetric matter, or for neutral 
current reactions this prescription results in the correct 
expression, but in asymmetric matter and for charged 
current reactions it is in error. To obtain the correct 
expression for the mean field polarization function from 
the free gas results of Q it is necessary to make both 
replacements given in Eq. 1161 



C. Correlations and Collisional Broadening 



FIG. 3: The top panel shows the total absorption inverse 
mean free path as a function of incoming neutrino energy 
for electron neutrinos (solid lines) and electron antineutri- 
nos (dashed lines). The dot-dashed line shows the effective 
bremsstrahlung inverse mean free path. In both panels the 
black lines include mean field effects and the red lines assume 
a free gas response function. The bottom panel shows the 
ratio of the total electron neutrino capture rate to the total 
electron antineutrino capture rate. Beta-equilibrium has been 
assumed and the temperature has been fixed at 8 MeV. 



approximation for electron antineutrino destruction. For 
electron antineutrinos at low energies, bremsstrahlung 
dominates the capture rate. Mean field effects push 
the energy region were bremsstrahlung is dominant to 
larger neutrino energies. This suggests that varying the 
assumed bremsstrahlung rate will also affect the spec- 
trum of the electron antineutrinos. In the bottom panel, 
the ratio of the electron antineutrino mean free path to 
the electron neutrino mean free path is shown as a func- 
tion of energy with and without the affect of mean fields. 
The large asymmetry induced between electron neutrino 
and antineutrino charged current interactions when mean 
fields are properly included is plainly visible. 

The formalism of [l| includes this effect, and was used 
to calculate the neutrino interaction rates employed in 
the models presented in [ll[ and in section IIIII of this 
work. However, the formulae in [26| and Q for charged 
current rates neglect the potential energy difference in 
the nucleon kinematics. In a procedure is advocated 
for including mean fields in which the effective chemi- 
cal potential, fit of each species is calculated from the 



In addition to the mean field energy shift, interactions 
correlate and scatter nucleons in the medium. The ex- 
citation of two or more nucleons by processes such as 
v e + n + n — > n+p + e~ and v e + n + p — > p + p + e~ alter 
the kinematics of the charged current reaction. Typically, 
these two-particle reactions introduce modest corrections 
to the single-particle response when the quasi-particle 
life-time is large. However, they can dominate when: (i) 
energy-momentum requirements are not fulfilled by the 
single particle reaction; (ii) final state Pauli blocking re- 
quires large energy and momentum transfer; (iii) or both. 
Such circumstances are encountered in neutron star cool- 
ing, where the reaction n — > e~ + p + D e is kinematically 
forbidden at the Fermi surface under extreme degeneracy 
unless the proton fractions x p > 10% [U[2^]. Instead, 
the two-particle reaction n + n — > e~ + p + D e , called 
the modified URCA reaction, is the main source of neu- 
trino production [29| . At temperatures encountered in 
PNS cooling, energy-momentum restrictions do not for- 
bid the single-particle interactions, but they do strongly 
frustrate them due to final state blocking. 

The excitation of two particle states in neutral cur- 
rent reactions has been included in a unified approach 
described in 0] and incorporated into the total response 
function by introducing a finite quasi-particle lifetime r. 
This naturally leads to collisional broadening allowing 
the response to access multi-particle kinematics and al- 
ters both the overall shape and magnitude of the response 
function [j| 0]- Here, as a first step, we adapt the gen- 
eral structure of the response function from [7| to show 
that two-particle excitations play an important role in the 
charged current process. We include a finite r through 
the following ansatz for the imaginary part of the polar- 
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which is obtained by replacing the energy delta-function 
in the Fermi gas particle-hole polarization function (see 
Eq. [4] and Eq. [6J by a Lorentzian with a width Y = 1 jr. 
Here, as before z = (qo+fi)/T and Ae p+q — t P + q ~ e p .The 
Lorentzian form is obtained in the relaxation time ap- 
proximation discussed in Q, and is valid when \qo\r 3> 1- 
The quasiparticle lifetime r is a function of the quasi- 
particle momentum, g, go and the ambient conditions. Its 
magnitude and functional form at long- wavelength is con- 
strained by conservation laws. For the vector-response, 
t — Y oo in the limit g — > due to vector current conser- 
vation. However, because spin is not conserved by strong 
tensor and spin-orbit interactions, the nucleon spin fluc- 
tuates even at q — > and the associated spin relaxation 
time T a is finite Since the spin response dominates 
the charged current reaction, in what follows we shall use 
Eq.[TH]only to modify the spin part of the charged current 
response. We, however, note that the multi-particle re- 
sponse in the vector channel warrants further study since 
the typical momentum transfer q ~ fi e is not negligible. 

For the spin relaxation time T a we use results calcu- 
lated in Ref. Q which indicate that it decreases rapidly 
with both density and temperature. The typical range 
of values of 1 /r ff obtained from [8j but including a 50% 
variation over their quoted values is shown in Figure [T] for 
conditions in the neutrino sphere region. Using these val- 
ues as a guide we study the effects of collisional broaden- 
ing on the v e and v e cross-sections. The differential cross- 
section for the axial portion of the process z/ e +n — > p+e~ 
is shown in Figure g] for T = 8 MeV, n B = 0.02 fm~ 3 , 
and Y e — 0.027. The initial neutrino energy is E Ve = 12 
MeV. As before the differential cross-section is plotted 
as function of the outgoing electron energy. The result 
with r — > recovers the single-particle response with the 
mean field energy shift included. Representative values 
of r = 1,2,4 MeV are chosen to approximately reflect 
the findings of Q for these ambient conditions. The col- 
lisional broadening seen in Figure 2] is quite significant. It 
increases the the axial portion of the cross-section by ap- 
proximately 20%, 44% and 80%, for T = 1,2,4 MeV, re- 
spectively. Together, the mean field energy shift and col- 
lisional broadening push strength to regions where elec- 
tron final state blocking is smaller resulting in an overall 
increase in the electron neutrino absorption rate. 

While mean field effects reduce the v e cross-section, 
collisional broadening will tend to increase it by accessing 
kinematics where —go is larger. This is shown in the 
inset of Figure 0] where the P e cross-section for the same 
ambient conditions and for £ , Pe =12 MeV is plotted as a 
function of the positron energy = E v — go . The units 
are arbitrary and the plots only serve to illustrate the 
relative effect of multi-pair excitations. We choose the 
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FIG. 4: The axial portion of the u £ (main panel) and v a (in- 
set) absorption cross-section including collisional broadening. 
This shifts a significant fraction of the response to larger E e 
where there is larger lepton phase space available. The ambi- 
ent conditions and neutrino energy are the same as those in 
figure The dashed lines show the RPA response, including 
both mean fields and collisional broadening. The dotted line 
in the inset panel shows the free gas response. 



same values of Y as for the v e case. Here broadening due 
to multi-pair excitations has a more significant effect than 
for v e absorption. However, despite this enhancement, 
the response that includes the mean field energy shift 
and collisional broadening is still much smaller than the 
free gas response. 

In addition to multi-pair processes, weak charge 
screening in the medium can also affect the charged cur- 
rent response. Screening due to correlations has been in- 
vestigated in the Random Phase Approximation (RPA), 
where specific long-range correlations are included by 
summing single-pair "bubble" or particle-hole diagrams. 
Additionally, this approach ensures consistency between 
the response functions and the underlying equation of 
state in the long wavelength limit. For charged currents, 
calculations reported in Q and [ij indicate that the sup- 
pression is density and temperature dependent. It can 
be as large as a factor of 2 at supra nuclear density, but 
at densities of relevance to the neutrino sphere where 
P ^ 10 13 g/ cm3 the corrections are « 20%. More im- 
portantly, the suppression found in 0] and [|| for the 
charged current rate is a weak function of reaction kine- 
matics and can viewed as a overall shift of the response 
in Fig. [5J aside from regions were significant strength is 
shifted to collective modes. The energy and momentum 
restrictions discussed previously apply also to the RPA 
response, and the mean field energy shift is important to 
include in the calculation of the particle-hole diagrams. 
They were included in [|| but omitted in Q . 

To include correlations between particle-hole (p-h) ex- 
citations due to residual interactions in the spin-isospin 
channel using the RPA, we employ a constant interac- 
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tion (independent of momentum and density) in the spin- 
independent and spin-dependent particle-hole channels 
given by V^ h = 2 fm~ 2 and = 1.1 fm -2 , respectively. 
The residual interaction in spin-independent channel is 
consistent with underlying equation of state. The resid- 
ual interaction in spin-dependent particle-hole channel is 
retrieved from analysis of the Gamow- Teller transition 
in finite nuclei [30( . The RPA response functions for this 
simple form of the p-h interaction are then given by 



S T {qo,q) 



1 - e- z 
2 

1 - e- 



Im 



n 



MF 



Im 



i - K pll n M F 
n r 



i - e h n r 



(19) 



(20) 



and the real and imaginary parts of the polarization func- 
tions satisfy the Kramers-Kronig relation, 



V f 

Re U(q ,q) = - — / du 



Im n(w, q) 



5o 



(21) 



RPA correlations also act to redistribute the strength 
of the response. The RPA response is shown in Figure 
131 The Gamow- Teller resonance is clearly visible in the 
curves that do not include large amounts of collisional 
broadening. The extent to which this affects the inverse 
mean free paths can be gauged from the results presented 
in Table |U 

While collisional broadening tends to increase both v e 
and v e cross-sections, RPA correlations decrease the v e 
cross-section and enhance the cross-section for v e . Given 
the simplicity of our model for the p-h interaction, these 
results only serve to capture the qualitative aspects of the 
role of correlations. They nonetheless demonstrate that 
changes expected are small compared to corrections aris- 
ing due to a proper treatment of mean field effects in the 
reaction kinematics. Hence, in the following discussion of 
PNS evolution and neutrino spectra, we set aside these ef- 
fects due to RPA correlations and collisional broadening, 
and calculate the neutrino interactions only including the 
mean field energy shifts calculated as described in |l|. 



III. PROTO-NEUTRON STAR EVOLUTION 

To illustrate the effect of the correct inclusion of mean 
field effects in charged current interaction rates, as well 
as the importance of the nuclear symmetry energy, five 
PNS cooling models are described here. The models have 
been evolved using the multi-group, multi-flavor, general 
relativistic variable Eddington factor code described in 
12] which follows the contraction and neutrino losses of 
a PNS over the first ~ 45 seconds of its life. These start 
from the same post core bounce model considered in [l2[ 
and follow densities down to about 10 9 gcm -3 . There- 
fore, they do not simulate the NDW itself but they do 
encompass the full neutrino decoupling region. 

One model was run using neutrino interaction rates 
that ignore the presence of mean fields, but are appro- 
priate to the local nucleon number densities (i.e. the 




4 6 
Time (s) 

FIG. 5: Top panel: First energy moment of the outgoing 
electron neutrino and antineutrino as a function of time in 
three PNS cooling simulations. The solid lines are the av- 
erage energies of the electron neutrinos and the dashed lines 
are for electron antineutrinos. The black lines correspond 
to a model which employed the GM3 equation of state, the 
red lines to a model which employed the IU-FSU equation of 
state, and the green lines to a model which ignored mean field 
effects on the neutrino opacities (but used the GM3 equation 
of state). Bottom panel: Predicted neutrino driven wind elec- 
tron fraction as a function of time for the three models shown 
in the top panel (solid lines), as well as two models with the 
bremsstrahlung rate reduced by a factor of four (dot-dashed 
lines). The colors are the same as in the top panel. 



re-normalized chemical potentials, were used but we 
set AU = 0). The equation of state used was GM3. 
This model was briefly presented in [l2[ • Another model 
was calculated that incorporated mean field effects in the 
neutrino interaction rates and used the GM3 equation of 
state. A third model was run using the IU-FSU equation 
of state and including mean field effects but with every- 
thing else the same as the GM3 model. Additionally, two 
similar models were run with the bremsstrahlung rates of 
[H reduced by a factor of 4 as suggested by [3l| . The neu- 
trino interaction rates in all five models were calculated 
using the relativistic polarization tensors given in [l[ with 
the weak magnetism corrections given in Q. 

In the top panel of Figure O the average electron neu- 
trino and antineutrino energies are shown as a func- 
tion of time for the three models with the standard 
bremsstrahlung rates. As was described in [12], includ- 
ing mean field effects in the charged current interaction 
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TABLE I: 1/A in m" 1 for matter in beta-equilibrium at T = 8 MeV and various densities and E Vc = E De — 12 MeV. The 
entries in the table follow the notation a (6) = a x 10 b . In the last two columns, T is considered to density dependent and the 
values used are taken from Figure [1] 



Density (fm 3 ) 


1/A (m" 


1 ): no MF 


mf (r 


= 0) RPA (r 


— 


0) MF (r > 0) RPA (r > 0) 


Tib = 0.020 


1/A„ e 


5.9 (-4) 


5.2 (- 


-3) 


2.1 (- 


-3) 


7.5 (-3) 


3.9 (-3) 




1/A* e 


3.5 (-4) 


2.7 (- 


-5) 


6.5 (- 


-») 


4.5 (-5) 


6.0 (-5) 


n B = 0.006 




7.7 (-4) 


1.6 (- 


-3) 


1.2 (- 


-3) 


1.8 (-3) 


1.3 (-3) 




1/A„ e 


2.4 (-4) 


1.4 (- 


-4) 


2.0 (- 


•4) 


1.4 (-4) 


1.5 (-4) 


n B = 0.002 


1/A,e 


5.3 (-4) 


6.5 (- 


-4) 


5.9 (- 


-4) 


6.8 (-4) 


6.1 (-4) 




l/Xv t 


1.5 (-4) 


1.3 (- 


-4) 


1.4 (- 


-4) 


1.3 (-4) 


1.3 (-4) 



rates significantly reduces the average electron neutrino 
energies because the decreased mean free paths (relative 
to the free gas case) cause the electron neutrinos to de- 
couple at a larger radius in the PNS and therefore at a 
lower temperature. Conversely, for the electron antineu- 
trinos the mean free path is increased, they decouple at 
a smaller radius and higher temperature, and therefore 
their average energies are larger. Mean field effects serve 
to shift the average neutrino energies by around 25% at 
later times. The antineutrino energies are also slightly 
larger than the values reported in [l2j because of the re- 
duced bremsstrahlung rate. 

To illustrate the properties of the region where neu- 
trino decoupling occurs, a snapshot of the decoupling re- 
gion as a function of neutrino energy is shown in Figure 
[51 In this work, the "decoupling region" is defined as 
the region where the Eddington factor /i = F g /N g obeys 
the condition 0.1 < A < 0.5. Here, F g is the neutrino 
number flux in energy group g divided by the speed of 
light and N g is the neutrino number density in energy 
group g (see [l2l]). This approximately defines the re- 
gion over which neutrinos transition from being diffusive 
to free-streaming. Higher energy electron neutrinos de- 
couple at a larger radius and therefore a lower density 
and temperature. At these radii, AU is smaller than the 
temperature and the inclusion of mean fields in the in- 
teraction rates should not significantly change the high 
energy electron neutrino mean free paths. At lower neu- 
trino energies, AU is significantly larger than the temper- 
ature in the decoupling region and the presence of mean 
fields strongly affects the opacity. As time progresses, 
the average neutrino energies become lower and decou- 
pling occurs in conditions at which mean field effects be- 
come increasingly important. Decoupling also occurs at 
a higher density for lower energy neutrinos, where both 
multi-particle processes and RPA corrections can poten- 
tially become important. 

Additionally, there are significant differences between 
the two models which include mean field effects but use 
different equations of state. As was described above, 
the GM3 equation of state has a smaller symmetry en- 
ergy than the IU-FSU equation of state at sub-nuclear 
densities and therefore has a smaller AU in the neu- 
trino decoupling region. This suggests that GM3 should 
have slightly larger electron neutrino average energies 




e v (MeV) 



FIG. 6: The thermodynamic conditions and nucleon potential 
difference characterizing the region where electron neutrinos 
decouple as a function of neutrino energy. These values are 
taken from the PNS model which employed the IU-FSU equa- 
tion of state at 3.3 seconds after core-bounce. At this point 
the average electron neutrino energy is 8.3 MeV. 



and slightly lower average electron antineutrino energies. 
The results of self-consistent PNS simulations are some- 
what more complicated than this simple picture, mainly 
because the equilibrium electron fraction near the neu- 
trino sphere also depends on the nuclear symmetry en- 
ergy which affects the charged current rates (see Figure 



9 



[T]). Still, there is a larger difference between the aver- 
age electron neutrino and antineutrino energies through- 
out the simulation (relative to GM3) when the IU-FSU 
equation of state is used, as expected. 

The moments of the escaping neutrino distribution 
along with the electron neutrino number luminosities can 
be used to calculate an approximate NDW electron frac- 
tion Hi 



Y, 



e,NDW 



N Ve (<7(e)„,0 



(22) 



where N are the neutrino number luminosities and (<r(e)) 
are the energy averaged charged current cross-sections in 
the wind region. The approximate NDW electron frac- 
tion as a function of time for the five models is shown in 
the bottom panel of Figure GO The low density charged 
current cross-sections from [33[ were used. First, it is 
clear from this plot that mean field effects significantly 
decrease the electron fraction in the wind. This is mainly 
due to the increased difference between the electron neu- 
trino and antineutrino average energies caused by the 
effective Q value induced by the mean field potentials. 
Second, increasing the sub-nuclear density symmetry en- 
ergy decreases the electron fraction in the wind. This 
in turn implies that nucleosynthesis in the NDW may 
depend on the nuclear symmetry energy because it is 
sensitive to electron fraction in the wind [e.g. 15]. Still, 
this effect is not particularly strong because the increase 
in the electron neutrino cross section for increased AU 
is partially mitigated by the larger equilibrium electron 
fraction predicted for models with a larger nuclear sym- 
metry energy. 



IV. CONCLUSIONS 

In this work, we have discussed the physics of charged 
current neutrino interactions in interacting nuclear mat- 
ter at densities and temperatures characteristic of the 
neutrino decoupling region in PNS cooling. Addition- 
ally, models of PNS cooling have been run to assess the 
importance of changes in the charged current rates to the 
properties of the emitted neutrinos. Our main findings 
are: 

• The mean-field shift of the nucleon energies al- 
ters the kinematics of the charged current reac- 
tions. Under neutron-rich conditions it increases 
the Q-value for v e absorption and decreases it 
for v e . Due to final state blocking (electron 
blocking for electron neutrino capture and neu- 
tron blocking for electron antineutrino capture), 
the increase in the Q value leads to an exponen- 
tial (exp (AU/T)) increase in the v e cross-section 
absorption and reduces the v e absorption cross- 
section by exp {—AU /T). 



• The formulae for the charged rates developed in Q 
and (2(| neglect these effects and the prescription 
for incorporating mean field energy shifts outlined 
in [H is inconsistent. 

• The nuclear symmetry energy at sub-nuclear den- 
sity plays a crucial role in determining the magni- 
tude of the difference between the mean field neu- 
tron and proton potential energies, and through its 
effect on the Q-values increases the difference be- 
tween the mean free paths of u e and v e . This sen- 
sitivity to the symmetry energy is potentially ex- 
citing since supernova neutrino detection and nu- 
cleosynthetic yields may be able to provide useful 
constraints. 

• Our 'preliminary work indicates that multi-pair ex- 
citations favor kinematics where final state electron 
blocking is small because the energy/momentum 
constraints present when only single particle-hole 
(p-h) excitations are considered are relaxed. This is 
analogous to the importance of the modified URC A 
process in neutron star cooling. In contrast to 
mean field effects, multi-pair excitations decrease 
the mean free paths of both electron neutrinos and 
electron antineutrinos. 

• Nuclear correlation effects treated in the RPA de- 
crease the z/ e cross-section and enhance the cross- 
section for z7 e . However, preliminary calculations 
using residual interactions consistent with equation 
of state or derived from Gamow- Teller transitions 
of finite nuclei suggest that the changes are much 
smaller than the proper inclusion of mean field ef- 
fects in the reaction kinematics. Although it is diffi- 
cult to determine from the limited and approximate 
calculations performed for this work, it seems most 
likely that multi-pair excitations and RPA correc- 
tions will bring the average electron neutrino and 
antineutrino energies somewhat closer to one an- 
other (relative to the case were only mean fields 
are included). 

• As was shown in 12], the changes to the charged 
current mean free paths induced by the correct in- 
clusion of mean fields decreases the average energy 
of the electron neutrinos and increases the average 
energy of the anti-electron neutrinos emitted dur- 
ing PNS cooling. The difference is relatively large, 
it significantly alters the predicted electron fraction 
in the NDW, and may have observable effects. This 
result has recently been independently confirmed 
by0. 

• We have also directly shown that increasing the 
value of the nuclear symmetry energy at sub- 
nuclear densities decreases the electron fraction in 
the neutrino driven wind. Therefore, NDW nucle- 
osynthesis may put some constraint on the poorly 
known density dependence of the nuclear symmetry 
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energy, or vice versa. This potential astrophysical 
constraint is in addition to those discussed in [35| . 
We emphasize that it may be hard to disentangle 
this from the effects of mult i- particle excitations, 
both on the charged current reactions themselves 
and on the (related) bremsstrahlung rate. This ef- 
fect is also partially compensated by the symmetry 
energy dependence of the beta-equilibrium electron 
fraction. 

• The reduced mean free path of v e is also likely to af- 
fect the de-lcptonization time of the proto-neutron 
star and may account for differences in time-scales 
observed in simulations performed using equations 
of state with different symmetry energies. 

Our work also shows that multi-particle excitations and 
correlations can alter the charged current response by as 
much as as factor of two at densities realized in the neu- 
trino decoupling region. However, our simple treatment 
has large uncertainty and warrants further study before 
we can make reliable predictions for the difference be- 
tween v e and v e spectra. Since this difference affects nu- 
cleosynthesis, collective neutrino oscillations, and is po- 



tentially observable from the high statistics expected for 
a galactic supernova neutrino burst, our study here iden- 
tifies that there is still much work to pursue both with 
respect to the charged current reactions and equation of 
state of neutron-rich matter in the neutrino decoupling 
region. 
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